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Method and device for the processing of interference in signals received by an 
array of sensors 
BACKGROUND OF THE INVENTION 
1. Field of the Invention 

The invention relates to a method and device for processing and 
eliminating the interference present in one or more signals received by a network of 
N sensors. 

It can be applied to the elimination of deliberate or involuntary 
interference occupying all or part of the spectrum of satellite signals received by GPS 
(global positioning system) receivers. 

The invention can be applied to improving interference processing 
methods in different signal-processing systems. 

It can also be used to get rid of deliberate or involuntary interference in 
signals received by standard receivers. 

2. Description of the Prior Art 

Systems for anti-interference processing in antenna networks presently 
use methods in which the entire band of the GPS signal received as input data is 
taken into account. 

In most of these methods, an apparent antenna is formed by the weighted 
combination of the signals coming from elementary sensors. What is done actually 
is to use a network of spatially separated sensors and, by a "constructive" or 
"destructive" combination, to highly attenuate the signal in all the directions 
identified as being occupied by one or more interference phenomena. Typically, 
these are standard principles of the CRPA (Controlled Radiation Pattern Antenna) 
implementing power inversion algorithms that are particularly well suited to useful 
noise signals whose level is lower than that of thermal noise, which is the case with 
GPS signals. 

To determine the coefficients of combination mentioned here above, the 
CRPA algorithm uses the principle described here below with reference to figure 1. 

The analog radioelectrical GPS signals, s ; , are received by the N sensors 
Ci of an antenna array. These signals Si have a spectrum constituted by a 20 MHz 
band centered on the frequency Li = 1.575 GHz (carrier frequency) and the 
frequency L 2 = 1.273 GHz, these two carrier frequencies being known in the GPS 
field. They are transmitted to a set 1 of transposition circuits, to be transposed to an 
intermediate frequency Fi lower than the carrier frequency Li (or L 2 ). The frequency 
transposition is achieved by methods known to those skilled in the art, such as the 
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methods described in the patent FR 2.742.612 by the present applicant for example. 
These signals thus taken to an intermediate frequency may be filtered. All the 
processing operations are performed by means of an analog process known to those 
skilled in the art. The filtered signals are then digitized by means of an ADC (an 
analog-digital converter) 2 that works at a chosen sampling frequency Fe to comply 
with the Shannon theorem. The ADC generates digital samples that contain GPS 
signals at a sampling frequency rate Fe and throughout the band of the useful signal, 
and are applied to a computation unit 3 and to a processing block 4. 

The computation unit 3 uses a CRPA type algorithm and a power 
inversion computation to identify the directions in which interference sources are 
present. This unit 3 determines the different weighting coefficients w ; to be applied 
to the digital samples. 

The weighting coefficients Wj are applied at input of the processing block 
4 to the samples Xi coming directly from the ADC 2, the unit 4 being adapted to 
making the interference sources disappear in the reconstituted samples, for example 
by combination of the weighted samples. 

The algorithm used to determine the weighting coefficients to be applied 
to the samples is especially well suited to signals known as narrowband signals, 
namely continuous wave (CW) type signals or signals with low frequency spread, 
typically for signals having a frequency^* to centerf reque ncy ratio that is smaller than 
unity. When interference comes into play on a wide frequency band, for example on 
the entire 20 MHz band in the case of the GPS P-code signal present in the L 2 band 
or again the C/A code present in the Li band, the interferences are less well 
eliminated by the power inversion algorithm. Or it is more likely that the number of 
degrees of freedom available, hence the number of interference phenomena that the 
receiver has to withstand, are thereby reduced. 

Furthermore, in the case of mobile carriers (GPS type receivers or 
stations comprising GPS receivers) and/or for mobile interferences in space, the 
estimation of the power and the combination to be made is more noise-affected. It is 
therefore less precise instantaneously and may result in phase leaps in the 
reconstituted GPS signal that will substantially disturb its nominal operation. In one 
of the methods used to overcome this problem, a smoothing stratagem is integrated 
into the processing algorithm. This is done, for example, by the addition of fictitious 
noise in order to reduce the resultant noise on the weighting coefficients and, 
therefore, on the phase of the resultant signal. Such stratagems may, however, 
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reduce the sensitivity of the anti-disturbance system, namely the level of minimum 
reference from which the power reversal algorithm will "perceive" and process the 
interference. The addition of fictitious noise raises the general floor above which the 
algorithm "perceives" the interference and above which the "minor" interference is 
not seen. 

SUMMARY OF THE INVENTION 

The object of the invention relates to a signal-processing method used to 
eliminate interference in a signal received by a network of N sensors, for example a 
satellite signal received by a GPS receiver. 

The object of the invention relates to a method to eliminate interference 
occupying at least one part of the spectrum of one or more signals received by a 
network of N sensors, the method comprising at least the following steps: 

• subdividing each sample Xj of signals into K frequency bands, 

• weighting the samples x ik obtained by subdivision, with weighting coefficients 
Wik determined by power inversion processing, 

• combining the different weighted coefficients w ik .Xi k by given frequency band 

N 

index k to obtain signals s k corresponding to ^w ik .x ik , and then carrying out the 

1=1 

combination of the signals s k for the totality of the bands K. 

The power inversion processing is, for example, of the CRPA type. 

The invention also relates to a method to eliminate interferences 
occupying a part of the spectrum of a signal received by a network comprising N 
sensors, wherein the method comprises at least the following steps : 

• digitizing the signals Sj received by the sensors in N digital samples x;, 

• transmitting the Xj digital samples to K filters G k in order to subdivide each 
sample xj into K frequency bands, 

• applying the Xj k samples obtained by subdivision to : 

• a computation unit adapted to determining the weighting coefficients w ik , by 
power inversion processing, 

• a processing block adapted to : 
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• combining the different weighted coefficients Wj k .x ik for a given filter 

N 

index k in order to obtain signals s k corresponding to ^ w tt jc a . 

• combining the signals s k in order to obtain a signal S' that is totally or 
mostly free of interferences. 

The object of the invention also relates to a device to eliminate the 
interferences in one or more signals Sj received by an array of N sensors comprising 
at least one set of means adapted to subdividing each sample Xi of signals into K 
frequency bands, weighting the samples x ik obtained by subdivision, combining the 
different weighted coefficients w ik .x ik by given frequency band index k to obtain 

N 

signals s k corresponding to ^w tk .x ik , combining the signals s k for the totality of the 
bands K. 

The power inversion processing is, for example, of the CRPA type. 
According to one embodiment, the device comprises at least: 

• one signal reception chain comprising circuits for the frequency transposition of 
the frequency of the initial signal to an intermediate frequency and an ADC to 
convert the signal S into N digitized samples, 

• a device adapted to subdividing each digitized signal x ; , into K frequency bands, 
in order to give N*K samples x ik , 

• a computation unit receiving the N*K samples and suited to determining 
weighting coefficients w ik , by power inversion processing, 

• a processing block receiving the weighting coefficients w ik and the samples x ik , 
said block being suited to the application of the weighting coefficients to the 
different samples, carrying out the combination firstly for a given index k of the 
x ik weighted samples with k of varying from 1 to K and secondly the K signals s k 
with k varying from 1 to K, in order to obtain a signal S'. 

The method and the device according to the invention are applied for 
example to eliminating the interferences in the signals sent by a satellite and received 
by a GPS receiver or again by a spread-spectrum positioning system or again a 
spread-spectrum navigation and communications system. 
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In particular, the invention offers the following advantages: 

• it very substantially strengthens the capacities of resistance to disturbing 
phenomena (deliberate or involuntary interference), 

• based on the principle of "network" processing to carry out "spatial" elimination, 
the invention is released from the need to carry out the classically used 
"narrowband" approximations, 

• it brings less noise into the narrowband for adaptive processing, thus tending to 
increase the sensitivity of the algorithm used, 

• by the addition of a Kalman filtering: 

• it absorbs the processing defects related to the dynamics of the carriers and the 
disturbing phenomena, and 

• gives an adaptive process for the correction of defects liable to be introduced 
by the hardware structure, the changes in the capacities of the components as a 
function of thermal phenomena for example. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the invention shall appear from the 
following the detailed description made with reference to the appended drawings, of 
which: 

• Figure 1 shows an exemplary prior art GPS receiver, 

• Figure 2 gives a diagrammatic view of the first implementation of the invention 

• Figure 3 shows a second implementation of the invention integrating a Kalman 
filter. 

MORE DETAILED DESCRIPTION 

In order to understand the object of the invention more clearly, the 
following description is given by way of an illustration that in no way restricts the 
scope of the invention for the processing of interference in signals received by GPS 
receivers. 

In a manner similar to that of figure 1, the device has an array of N 
sensors Ci, a frequency transposition block and an analog-digital converter not 
shown in figure 2 for reasons of simplification. 

The N samples coming from the ADC 2 (figure 1) are applied to a device 
20 adapted to carrying out a frequency subdivision. The frequency subdivision is 
performed by using a set of K finite impulse response (FIR) digital filters. The 
device 20 is provided with N input channels 20; corresponding to the N samples Xj, i 
being an index designating a sample, and N*K output channels 20 ik , with k being the 
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index corresponding to the filter used. A sample Xi is applied to the K filters Gk so as 
to obtain K the digital signals designated by Xik, corresponding to K bands narrower 
than the initial band of the signal. 

The characteristics of each and/or of all of the K filters Gi to Gk are 
chosen so that the sum of the frequency band thus obtained for each sample x,k, 
reconstitutes the total useful band fully or as fully as possible. Each sample has a 20 
MHz useful band corresponding to the useful band of the GPS signal received on the 
sensor Ci. 

The band separation process is achieved preferably by digital means. 
This provides for a precise adjustment of the coefficients of the different filters in 
order to obtain distortion-free reconstitution of the total band. 

The samples xjk thus obtained are applied firstly to a computation unit 21 
and secondly to a processing block 22. 

The computation unit 21 is programmed to carry out a CRPA type power 
inversion processing and compute the dedicated weighting coefficients Wik, band by 
band for the N*K samples. At the end of this computation, the method is in 
possession of K sets of weighting coefficients (N*K coefficients), to be applied to 
the different samples Xik for example at input of the processing block 22. The 
weighting coefficients thus obtained are better suited to the elimination of the K 
potential interference bands. 

The processing block 22 is adapted to combining the weighted samples 
Wik.Xik. The combination step is carried out for example by initially combining the 
different weighted samples for a given filter index k, in achieving a variation in index 

N 

i from 1 to N, to obtain several signals Sk corresponding to ^jw jk .x ik . In a second 

K 

stage of this combination step, the signals Sk, ^ s k are summed. The sum represents 

the reconstituted signal S' exempt or practically exempt from interference. The 
different computations are made by means of appropriate processing algorithms, and 
the components used could be of the FPGA or ASIC type. 

Advantageously, this embodiment is used to overcome the "narrowband" 
limitation of the commonly used CRPA type adaptive methods of power inversion. 
Furthermore, by working on narrower bands then the initial signal band, the noise 
level is reduced to the processing level. Hence, for equivalent filtering processing, 
the sensitivity of the method is increased. 
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Figure 3 describes a second exemplary embodiment of the invention 
where the similar elements taken up in figure 2 relate to the same references. This 
embodiment is especially well suited in the case of mobile interference or mobile 
carriers. 

In this example, the N*K weighting coefficients obtained by the power 
inversion computation are applied in a dynamic filtering step, by using for example a 
Kalman type filter 30. The filter made by means of an adapted device, has the 
function especially of separating the directional coefficients from the N*K 
coefficients (with a high dynamic range or related to the dynamic range of the 
disturbing phenomena) from the distortions related to the reception lines (continuous 
components on a distant horizon). 

The dynamic range of the disturbing phenomenon is, for example, 
spectral, of the spectral sweep jamming type or again it may be an loaded type of 
geographical jammer. Again it may be disturbance from jammer switching or it may 
be a pulsed jammer type of temporal disturbance. 

By adapting the Kalman filter to the different dynamic ranges, it is 
possible to resorb a part of the problems of dynamic range related to the tracking of 
interference during a movement, for example a severe operational constraint while, at 
the same time, correcting the receiver distortions, such as HF defects in particular: 
phase matching, amplitude etc, which limit the performance of the elimination. 

Classically, in a Kalman filter, the matching is done by the judicious 
choice of the « model noise ». This noise is generally fixed and is defined at the 
designing stage but may also be defined as a function of criteria that do not arise out 
of the measurements found. 

The filtered coefficients are then sent to the processing block 22 to 
combine the different weighted samples. This operation is carried out by frequency 
band, as described with reference to figure 2. 

The total signal after processing is then reconstituted by summation, for 
example before it is used according to the known prior art methods as a signal 
obtained by a standard CRPA operation. 

Without departing from the framework of the invention, the method can 
be applied in the field of inertia/GPS hybridization and also in any field used to 
separate the dynamic values included in the weighted coefficients. 

The method can also be applied to all the signals of a spread-spectrum 
positioning system such as the GPS, the GLONASS (Global Orbiting Navigation 
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Satellite System), Galileo or any other spread-spectrum navigation and 
communications system. 



